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Abstract 

A search for Z)‘^/o(980) decays is performed using 3.0 fb“^ of pp collision 

data recorded by the LHCb experiment during 2011 and 2012. The /o(980) meson 
is reconstructed through its decay to the tt+tt" final state in the mass window 
900MeV/c^ < m(7r+7r“) < 1080MeV/c^. No significant signal is observed. The first 
upper limits on the branching fraction of B{Bg —>■ 71°/o(980)) < 3.1 (3.4) x 10“® are 
set at 90 % (95 %) confidence level. 
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1 Introduction 


Understanding the qnark-level snbstrnctnre of the scalar mesons is one of the main 
challenges in hadronic physics. The nnmber of observed states, and their masses and 
branching fractions, snggest that there is a contribntion from fonr-qnark wavefnnctions 
in addition to qq, and possible gluonic, degrees of freedom [^[^. However, the extent of 
mixing between the different components is nnclear. 

Measurement of the relative production of scalar mesons in and 5° meson decays can 
help to address this issue Measurements of 5°^^ —)■ J/V’/ decays, where / represents 

either the /o(500) (also known as a) or the /o(980) meson, and / —)■ 7 r'*~ 7 r ~ |5| -[8] have 
already provided important insight into the structure of the scalar mesons (^|10]T Studies 
of D^f decays provide complementary information to the 5°^^ —)■ J/V’/ case Jll| . 

Measurements of the branching fractions of B^ —>■ D^f q{500) and B^ —)■ D^ fo{9 80) decays 
have been obtained from Dalitz plot analyses of B^ —)■ D^7r~^7r~ decays @1 [l^ , but there 
is no experimental result to date on the decays. 

In addition, under the assumption that the /o(980) meson has a predominant ss 
component, the B^ —)■ D^fo{980) decay mode can be used to determine the angle 7 of 
the CKM unitarity tri angle fl4|[I^ , using the same methods that are applicable for the 
Bg —)■ D^(j) decay mode 
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Since the B^ —)■ D^(j) decay has recently been observed 
a signal for the B^ —)• D^f q(980) channel is expected if the branching fractions of the 


two decays are comparab 

X 10-5 “ 
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e. 


An explicit calculation predicts B Zi)°/o(980)) = 


In this paper, the result of a search for the B^ —)■ D^f q{980) decay is presented. The 
inclusion of charge conjugated processes is implied throughout the paper. The hnal 
state is reconstructed through the —)■ and /o(980) —)■ tt+tt” decays. The 

decay-time-integrated branching fraction is measured under the assumption that the 
Bg —)■ decay proceeds uniquely via the /o(980) resonance within the selected 

mass window, 900MeV/c^ < m{7r~^7i~) < 1080MeV/c^. This approach was used for the 
hrst observation of B^ —)■ J/'0/o(98O) decays it is also justihed by the fact that 


no other contribution to B^ —)• D^7r~^7i decays, for example through the B^^D 


TT 


process 24 , is expected at the current level of sensitivity. A further assumption is that 


the contribution from 5° — D^fo{980) decays, which is suppressed by the ratio of CKM 
matrix elements \V^fjV*g/ ~ 0.1, is negligible. Formally, the measurement is 
of the decay-time-integrated sum of the branching fractions for B^ — fo(980) and 
5° ^ T)°/o(980) decays. 

The analysis is based on 3.0 fb~^ of LHC pp collision data collected with the LHCb 
detector, with approximately one third taken at a centre-of-mass energy of 7TeV during 
2011 and the remainder at 8 TeV during 2012. The measurement is obtained by evaluating 
the ratio of branching fractions 


^ Zl°/o(980)) 


N{B^^ T)0/o( 980)) e(H0 ^ 

X 


_ X — flj 

6(50 ^ 50/o(980)) // 


5 ( 5 ° ^ DH+n-) N{B^ DH+n-) 
from which the absolute branching fraction for 5° —)■ D^fo(980) decays is determined 
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using the known value of B{B^ —)■ 13 . The yields N{Bg —)■ Z)°/o(980)) and 

N{B^ —)■ are determined from separate extended maximnm likelihood hts to 

the distributions of selected D^fo{980) and D^7i~^7i~ candidates in both the B candidate 
mass and the output of a neural network (NN) used to separate signal from combinatorial 
backgronnd. The combined reconstruction and selection efficiencies, e(i?° —?■ Zi)°/o(980)) 
and e(-B° —?■ are determined from simulated samples with data-driven corrections 

applied. The ratio of fragmentation fractions inside the LHCb acceptance has been 
measured to be fs/fd = 0.259 ± 0.015 (^. Equation 0 corresponds to a branching 
fraction for /o(980) —)■ tt+tt” of 100%, which is the conventional way to quote results for 
decays involving /o(980) mesons. 


2 LHCb detector 


The LHCb detector |26[|27| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing 6 or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 28 


surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream of the magnet. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of 
a track to a primary vertex, the impact parameter (IP), is measured with a resolution of 
(15 -|- 29/pt) M-ni, where px is the component of the momentum transverse to the beam, 
in GeV/c. Different types of charged hadrons are distingnished using information from 
two ring-imaging Gherenkov detectors 30 . Photon, electron and hadron candidates are 


identihed by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Mnons are identihed by a 
system composed of alternating layers of iron and multiwire proportional chambers 


31 


The trigger 32 consists of a hardware stage, based on information from the calorimeter 


and muon systems, followed by a software stage, in which all tracks with px > 500 (300) MeV 
are reconstructed for 2011 (2012) data. The software trigger requires a two-, three- or 
four-track secondary vertex with signihcant displacement from the primary pp interaction 
vertices (PVs). At least one charged particle must have px > 1.7 GeV/c and be inconsistent 
with originating from a PV. A multivariate algorithm 33 is used for the identihcation of 


secondary vertices consistent with the decay of a 6 hadron. 

Simulated events are used to characterise the detector response to signal and certain 


types of background events. In the simulation, pp collisions are generated using Pythia 34 
with a specihc LHGb conhguration 


35 


Decays of hadronic particles are described 

The 


in which hnal state radiation is generated using Photos 37 


by EvtGen 

interaction of the generated particles with the detector and its response are implemented 
using the Geant 4 toolkit 


38 as described in Ref. 39 
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3 Selection 


Candidates consistent with the decay chain —)■ D^7i~^7r~ with —)■ K^7t~ are 

selected. The selection procednre involves applying a preselection to the data sample 
before using a NN to reduce the combinatorial background. The NN is trained with the 
preselected data sample, using the sPlot method 41 with the B candidate mass 


as discriminating variable to separate statistically the signal and background categories. 
The input variables to the NN are related to the kinematic properties of the candidate, 
its isolation from the rest of the pp collision event, and the topology of the signal decay 
chain. Full details of the preselection and NN training can be found in Ref. 42 . The 


four hnal-state tracks must also satisfy particle identihcation (PID) requirements. Signal 
candidates are retained for further analysis if they have invariant mass in the range 
5100-5900 MeV/c^. A requirement that the NN output is greater than —0.7 removes 77% 
of combinatorial background and retains 95 % of —)■ Zl°/o(980) decays. 

A requirement m{D^7T~) > 2.10 GeV/c^ is used to remove candidates that predominantly 
originate from B^ — D*~n~^ decays. A further requirement, m{D^7i~^) < 5.14GeV/c^, 
is used to remove backgrounds from B~^ —)■ decays combined with a random vr” 

candidate. This source of combinatorial background is kinematically excluded from the 
signal region, but causes structure in the mass distribution at higher B candidate mass. A 
similar contribution from R+ —decays cannot be vetoed in the same way, and 
must therefore be considered further as a source of background. 

Following all selection requirements, approximately 1 % of events contain more than 
one candidate. All candidates are retained for the subsequent analysis; the associated 
systematic uncertainty is negligible. 


4 Determination of signal yield 

The yields of B^ —>■ Zl°/o(980) and B^ —)■ D^7T~^7r~ decays are obtained from two separate 
extended maximum likelihood hts to the distributions of NN output and B candidate mass 
for selected candidates. The only difference between the samples used in the two hts is that 
the former has an additional requirement of 900MeV/c^ < m(7r^7r“) < 1080 MeV/c^. The 
yield of B^ —)■ decays in the latter ht is expected to be negligible compared to 

the large yields of B^ decays and combinatorial background, and is therefore hxed to zero. 
However, a signihcant number of B^ —)■ D^7r~^7T~ decays are expected to remain within the 


/o(980) mass window 13 , and therefore both B^ and B^ components are included in the 
former ht. 

The data are divided into hve bins of the NN output variable, dehned as [—0.70, 0.03], 
[0.03, 0.54], [0.54, 0.77], [0.77, 0.88] and [0.88, 1.00], and referred to hereafter as bins 1 to 
5, respectively. The hve bins contain a similar proportion of signal decays and increase in 
purity from bin 1 to bin 5. This choice of binning has been found to enhance the sensitivity 
whilst giving stable ht performance. 

The hts include components due to signal and combinatorial background as well as 
from partially reconstructed and misidentihed 6-hadron decays. The signal invariant mass 
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distribution is described by the sum of two Crystal Ball (CB) functions, with a shared 
mean and tails on opposite sides described by parameters that are hxed to values found in 
fits to simulated samples. 

The combinatorial background is modelled with the sum of two components. The hrst 
has an exponential shape, described by a parameter that is the same in all NN bins. The 
second originates from decays combined with a random pion candidate, and 

is modelled using a non-parametric shape determined from simulation. The limited sizes 
of the simulated samples used to obtain this and similar shapes are sources of systematic 
uncertainty. 


Partially reconstructed backgrounds occur from —)■ decays, with D*^ — 

and D*^ —> where the neutral pion or photon is not associated with the 

candidate, and from B^ —)■ decays where one tt’*' is also not associated with 

the candidate. The invariant mass shapes of these backgrounds are described with non- 
parametric functions derived from simulation. A global offset of the shape of the partially 
reconstructed background is determined from the ht to data to allow for differences between 
data and simulation 44 . 

Backgrounds from misidentihed 6-hadron decays arise from B^ —)■ and 

BO b{*)ok- 'Ti^ (hereafter collectively referred to as —?■ decays where 

the kaon is misidentihed as a pion and from —)■ decays where the proton is 

misidentihed as a pion. Simulation is used to obtain non-parametric descriptions of the 
invariant mass shapes. To obtain these shapes, the latest knowledge of the phase-space 
distributions of the decays (4^[4^[4^ , of the relative branching fractions of the B^ and 
B^ —)■ D^Kn modes J45| , and of the relative branching fractions of the decays involving 
and D*^ mesons (sT^is used. Data-driven estimates of the misidentihcation probability 
as a function of particle kinematic properties are also included. The relative yields in the 
NN bins are taken to be the same as for the signal decays. 

A total of 25 parameters are determined from the ht to the sample. These 

include yields of B^ —)■ D^tt^ 


TT 


partially reconstructed background, and the B^^^ 


decays, the total combinatorial background, the total 
—)■ and A° — D^*^^p7r~ misiden¬ 
tihed backgrounds. For B^ —)■ decays, combinatorial and partially reconstructed 

backgrounds, the fractional yields ft of each component in bins 1-4 are also free parameters, 
with the fraction in bin 5 determined as /s = 1 — fi- addition, the exponential slope 
parameter of the combinatorial background, the fraction of the combinatorial background 
from B~^ — D*^7r~^ decays, the fraction of the partially reconstructed background from 


B° ^ D 


*0 


vr^TT 


decays and the offset parameter of the partially reconstructed background 


are determined by the ht. Parameters of the signal invariant mass shape (the peak position, 
the width of the core CB function, and the relative normalisation and ratio of the CB 
widths) are also allowed to vary. Results of this ht are shown in Fig. 

The ht to the Zl°/o(980) subsample includes the same components as the B^ — D^n~^7r~ 
ht, with the addition of a second signal component to account for the possible presence of 
both B^ and B^ decays. The mass difference between the B^ and B^ mesons is hxed to the 
known value [^. The shapes for the B^ and B^ components are otherwise identical in both 
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Figure 1: Invariant mass distribution of candidates in the data sample with fit results 

overlaid, shown on a logarithmic scale. The components are as detailed in the legend. The labels 
(a) to (e) show the NN bins with increasing pnrity. The NN binning scheme is described in Sec. 01 


invariant mass and NN output. The following parameters are fixed to the values found in 
the —)■ ht: the fractional yields /* for the signal and partially reconstructed 

background components; the relative normalisation of the two CB functions; the ratio 
of widths of the CB functions; the fraction of the partially reconstructed background 
from decays and the offset parameter of the partially reconstructed 

background. In addition, the relative yields of the misidentihed background components 
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Table 1: Yields from the fit to the and Z)'^/o(980) samples. 




D^M980) 

R° ^ D^n+n- 

42 636 ± 362 

3 998 ± 87 

RO ^ R°/o(980) 

— 

29 ± 17 

Combinatorial 

90 150 ± 481 

11 064 ± 145 

Partially reconstructed 

50 950 ± 493 

3 759 ± 88 

R°) ^ DAAKtt 
yjo ^ dA)^ptt- 

9 225 ± 504 

852 ± 128 

4 923 ± 415 

154 ± 135 


The remaining 14 parameters are: the yields for —)■ D^fo{980) decays, D^7i~^7r 

decays, combinatorial and partially reconstructed backgrounds and for the B^^-^ —>■ D^*^^K7i 
and —>■ pn misidentihed backgrounds; the fractional yields of the combinatorial 

background in NN output bins, the exponential slope parameter of the combinatorial 
background and the fraction of the combinatorial background from B^ —decays; 
and the signal peak position and core width. Results of this £t are shown in Fig. 

The yields from the hts to the Zl°/o(980) and data samples are summarised in 

Table In total, 29 ± 17 5° —Zl°/o(980) decays are found, with a statistical signihcance 
of 2.2(7 obtained from \/—2 A In £, where Aln£ is the change in log likelihood from the 
value obtained in a £t with zero signal yield. 

5 Systematic uncertainties 

The systematic uncertainties that affect the ratio of branching fractions are summarised in 
Table Various effects contribute to the systematic uncertainties on the invariant mass 
£t and efficiencies, as described below. 

The tail parameters of the signal components for —>■ D^f q{980) and B^ —)■ 

decays are varied within the uncertainties from the ht to simulated events. For the 
B^ — )■ Zl°/o(980) ht, the relative normalisation and ratio of widths of the CB functions are 
varied according to the uncertainties from the ht to the B^ —mode. Combined 
in quadrature these contribute 11.5% to the systematic uncertainty. The systematic 
uncertainty from assuming that the NN response is identical for B^ —)■ Zl°/o(980) and 
B^ —)■ D^7T^7r~ decays is evaluated by correcting the fractional yields found in the 
B^ —!■ ht by the ratio of fractional yields found in simulated samples. This 

contributes 0.3 % to the systematic uncertainty. 

A second-order polynomial function is used to replace the exponential shape for the 
combinatorial background in both hts, giving a systematic uncertainty of 8.4%. Varying 
the smoothing of the non-parametric shape for B~^ —>■ decays gives the largest 

source of systematic uncertainty of 23.1 %; the size of this ehect is determined by that of 
the simulated background sample. The fractional yields of the B~^ — D*^7r~^ component 


from RO ^ DAAK+71- and R° ^ D^^k- '7r~^ decays are hxed to the expected value 


45 
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m{D n^Tt ) (MeV/c^) 


Figure 2: Invariant mass distribution of candidates in the -D*^/o(980) data sample with fit results 
overlaid, shown on a logarithmic scale. The components are as detailed in the legend. The labels 
(a) to (e) show the NN bins with increasing purity. The NN binning scheme is described in Sec. 01 


of the combinatorial background are fixed to the values found in simulation, rather than 
using the same fractional yields as the rest of the combinatorial background. This leads to 
a systematic uncertainty of 1.0 %. 

The smoothing of the non-parametric functions for —)■ T)*°7r+7r“ and —)■ 

is varied in both hts. Additionally, in the B'l — D^fo{980) sample, the relative 
normalisation of the shapes is varied within uncertainties from the value found in the 
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Table 2: Summary of systematic uncertainties on the ratio of branching fractions. 


Source Value 

Signal shapes 11.5% 

Combinatorial backgronnd shapes 24.6% 

Partially reconstrncted background shapes 6.9% 
Misidentihed background shapes 10.8% 

Efficiencies 2.5% 

Fragmentation fraction fs/fd 5.8% 

Total 30.7% 


—)■ D^'n'~^7T~ ht. Combined in quadrature these contribute 6.2 % to the systematic 
uncertainty. Allowing the fractional yields of the partially reconstructed background 
to vary in the Z)°/o(980) ht, instead of being hxed to values found in the ht, 

contributes 3.1 % to the systematic uncertainty. 

The misidentihed background shapes are also varied by changing the smoothing applied 
to the non-parametric function. Additionally, the simulation is not reweighted to the 
known phase-space distributions and the relative normalisation of the 5°^^ —?• 
shapes is varied within uncertainties. Together these contribute 6.7% to the systematic 
uncertainty. Corrections, derived from simulation, are applied to the fractional yields 
for the misidentihed backgrounds, which are assumed to behave like signal decays in the 
default ht. The sum in quadrature of the individual contributions gives a systematic 
uncertainty of 8.5 %. 

Potential biases in the ht procedure are investigated using an ensemble of pseudoexper¬ 
iments. Each of the pseudoexperiments is htted with the same ht model used to describe 
the data samples. This study shows that the ht is stable and well behaved and that the 
associated systematic uncertainty is negligible. 

The uncertainty on the ratio of reconstruction and selection efficiencies for the —)• 
-D°/o(980) and B^ —)■ hnal states contributes 2.5% to the systematic uncertainty. 

This includes statistical uncertainty from the sizes of the simulated samples as well as 
effects related to the choice of binning in kinematic variables in the evaluation of the PID 
efficiency and potential differences in the response of the hardware trigger. The simulated 
sample of B^ — i7°/o(980) decays is generated using a relativistic Breit-Wigner function 
with a width of 70 MeV for the /o(980) meson. The true lineshape of the /o(980) meson can 
differ from the assumed shape in a process-dependent way, which can affect the fraction 
of /o(980) —)■ 7r+7r“ decays that fall inside the selected m(7r+7r“) window. No systematic 
uncertainty is assigned due to this choice of /o(980) lineshape. Other possible sources of 
uncertainty on the ratio of efficiencies are negligible. 

The limited knowledge of the ratio of fragmentation fractions, fs/fd = 0.259±0.015 25 , 
contributes 5.8% to the systematic uncertainty. Combining all of the above sources in 
quadrature, the total systematic uncertainty on the ratio of branching fractions is found 
to be 30.7%. 
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6 Results and summary 


The relative branching fraction of —)■ D^fo{980) and —)■ D^7r^n~ decays is de¬ 
termined by correcting the ratio of yields for the relative efficiencies and fragmentation 
fractions, as shown in Eq. The total efficiencies are found to be e(i?° —)■ T)°/o(980)) = 
(0.76±0.02) % and e(i?° —)■ = (0.57±0.02) %. These values include contributions 

from the LHCb detector acceptance and from selection, trigger and PID requirements. 
The selection and trigger efficiencies are calculated from simulated samples with data- 
driven corrections applied. The PID efficiency is measured using a control sample of 
D*~ D^7r~, —)■ K^'k~ decays. Variation of the B^ —efficiency over the 

Dalitz plot is taken into account by weighting the simulation according to the observed 


Dalitz plot distribution 13 

Using Eq. Q the ratio of branching fractions is determined to be 


B{B^, ^ Dyo{9SO)) 

B{B0 D%+7r-) 


= (2.0 ±1.1 ±0.6) X 10 


-3 


where the hrst uncertainty is statistical and the second systematic. This result is obtained 
under the assumption that the B^ —)■ D° 7 r+ 7 r“ decays proceed uniquely via the /o(980) 
resonance within the range 900MeV/c^ < < 1080MeV/c^; no systematic uncer¬ 

tainty is assigned due to this assumption. The result can be converted into an absolute 
branching fraction by multiplying by B{B^ = (8.46 ±0.51) x 10“'^ 


13 


to give 


B{B° T)°/o(980)) = (1.7 ± 1.0 ± 0.5 ± 0.1) x 10 


-6 


where the third uncertainty is from B{B^ 
signfficant, upper limits of 


—)■ D^7r~^7i 


Since the signal yield is not 


B{B^^ T)°/o(980)) < 3.1 (3.4) x 10 


-6 


are set at 90 % (95 %) conhdence level. The statistical likelihood curve obtained from the 
fit is convolved with a Gaussian function of width equal to the systematic uncertainty. The 
limits obtained are the values within which 90 % (95 %) of the integral of the likelihood in 
the physical region of non-negative branching fraction are contained. 

In summary, a search for the B^ —)■ D^fo{980) decay has been performed using 3.0 fb~^ 
of pp collision data recorded by the LHCb detector in 2011 and 2012. No signfficant 
signal is observed, and a limit is set on the branching fraction that is below the predicted 
value . The small yield suggests that much larger data samples will be necessary in 
order to determine the angle 7 of the CKM unitarity triangle with B^ — Zl°/o(980) decays. 
Table ^ shows the current experimental status of measurements of the B^^^ —)■ Zl°/o(500) 
and N —)■ D^f q(980) branching fractions. The pattern of branching fractions is very 


(«) 

different to that for the B^^^ 


—)■ J/'ijj fo{500) and B^^-^ —)■ J/0/o(98O) modes 


results may provide insight into the substructure of the scalar mesons. 


These 
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Table 3: Results for branching fractions for T>°/o(500) and ^ Z)°/o(980) decays. 

All quoted results correspond to branching fraction for /o —?■ 7r'*'7r“ of 100%. There is no 
experimental result for B{B^ —^ Zl°/o(500)). 


/o(500) 

/o(980) 


g(go ^ D°/o)]T3} 


(11.2 ± 0.8 ± 0.5 ± 2.1 ± 0.5) X 10-^ 

(1.34 ±0.25 ±0.10 ±0.46 ±0.06) 


1 n-5 


^ D^fo) 


(1.7 ±1.0 ±0.5 ±0.1) X 10-6 
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